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Abstract. 

Luminous Compact Galaxies are enigmatic sources by many aspects. They can reach the luminosity of the Milky Way within 
a radius of only a few kpc. They also represent one of the most rapidly evolving populations of galaxies since they represent 
up to 1/5 of the luminous galaxies at redshift z = 0.7 while being almost absent in the local Universe. The measurement of 
their dynamics is crucial to our understanding of LCGs since this has the potential of telling us which physical process(es) 
that drives them, and ultimately to link them to the existing present-day galaxies. Here we derive the 3 dimensional velocity 
fields and velocity dispersion (cr) maps of 17 Luminous Compact Galaxies selected from the Canada France Redshift Survey 
and the Hubble Deep Field South with redshifts ranging from z=0.4 to z=0.75. We find that only 18% of them show rotational 
velocity fields typical of rotating disks, the others showing more complex kinematics. Assuming that LCGs are not too far from 
equilibrium, about half of LCGs then appear to be either non-relaxed objects, or objects that are not supported by velocity 
dispersion alone. This supports the view that an important fraction of LCGs are probably mergers. It brings additional support 
to the "spiral rebuilding scenario" in which LCGs correspond to a previous or post-merger phase before the disk re-building. 

Key words. Galaxies: evolution. Galaxies: formation , Galaxies: kinematics and dynamics , Galaxies: bulges 



1. Introduction 

Luminous Compact Galaxies (LCGs) (Mb < -20, ri,aif < 5 
/ijQ kpc and EWo{[OII]) >15A) correspond to the most rapidly 
evolving population seen in the UV (Lilly et al. 1998): they rep- 
resent ~ 20% of the 1 > z > 0.4 galaxies (Zheng et al. 2005), 
~50% of the emission line galaxies (see Paper I, Flores et al. 
2006) and almost vanish in the local Universe, with their num- 
ber density decreasing by factors 7-10 (Jangren et al. 2004; 
Garland et al. 2003; Werk et al. 2004). Moreover, LCGs con- 
tribute to 40-50% of the increase in the cosmic Star Formation 
Rate (SFR) density between z-0 and 1 as measured from rest 
frame UV luminosities (Lilly et al. 1995; Guzman et al. 1997) 
and 25% in the IR (Zheng et al. 2004). The spectra of LCGs 
reveal a mixture of old, intermediate and young stellar popu- 
lations (Hammer et al. 2001). Apart from their compactness, 
LCGs have properties surprisingly similar to those of other, 
more extended luminous IR galaxies and starbursts: They show 
a similar mix of stellar populations, extinction property distri- 
butions, stellar masses and SFRs (Hammer et al. 2005). This 
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led Hammer et al. (2001) to propose that LCGs are the progen- 
itors of present-day bulges of early type spirals. 

Recently, Hammer et al. (2005) proposed the so-called 
"spiral rebuilding scenario" to explain the formation of the 75% 
most massive local spirals. This scenario is composed of 3 ma- 
jor phases: a "pre-merger phase" during which two distant spi- 
rals merge, the "LCG phase" where all material from the pro- 
genitors fall into the mass barycenter of the system and form 
a bulge, and the "disk growing phase" where subsequently ac- 
creted material form a rotating disk. This scenario is partly sup- 
ported by Ostlin et al. (1999 and 2001) who obtained very com- 
plex Ha velocity fields for local LCGs, characteristic of what 
is expected from merging galaxies. Another alternative is the 
one proposed by Barton and Van Zee (2001): comparing HI 
and optical emission line widths of nearby LCGs candidates 
they argued that interactions and minor mergers of disk galax- 
ies may cause apparently compact morphology leading too to 
the formation of a bulge. Both views are not incompatible and 
can occur during different stages of the "LCG phase" described 
by Hammer et al. (2005). 

Compact galaxies have been extensively studied in the past. 
Koo et al. (1995) were the first to propose that some compact 
galaxies could be the progenitors of local dwarf ellipticals (dE), 
assuming that they experience a dramatic event of star forma- 
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tion before fading away by up to 5 magnitudes. Guzman et 
al. (1997) established a distinction between two types of com- 
pact galaxies: 60% present properties characteristic of young 
star-forming HII galaxies (in e.g. velocity widths, SFRs and 
mass-to-light ratios), whereas the remaining 40% constitute a 
more heterogeneous class of evolved starbursts, similar to lo- 
cal starburst disk galaxies. Philips et al. (1997) then suggested 
that the Hll-like compact galaxies are the best candidates to 
evolve into dEs. In this paper we study the most Luminous 
and most actively star forming fraction of the population of 
Compact Galaxies (LCGs), i.e. those which contribute most to 
the increase in the star formation rate density. It is important 
here to stress that the sample presented in this paper corre- 
spond to the brightest 25% of galaxies in the sample of Blue 
Compact Galaxies (BCGs) in the Hubble Deep Field (HDP) at 
0.1 < z < 1.3 studied by Guzman et al. (1997)', and would 
have mostly been classified as SB-disk like compact galaxies 
following their criteria. They thus do not correspond to the 
sub-class of compact galaxies for which Guzman et al. (1997) 
would have considered as possible progenitors of dEs. Notice 
that for another sub-class of compact galaxies. Philips et al. 
(1997) found that "one possibility is that they are disks form- 
ing from the center outward, and so the radius of the luminous 
material and enclosed mass are small compared to present-day 
spirals". 

Whereas the LCGs are important for understanding galaxy 
evolution since z = 1, their nature is thus still enigmatic. The 
aim of this paper is to evaluate the nature of their kinematics 
from a survey of 17 LCGs randomly selected from the CFRS 
and the HDF-South (HDFS) field at redshifts from z =0.4 to 
0.75, and to investigate if their dynamics is supported either by 
rotation or velocity dispersion. This will help in distinguishing 
a merger scenario from a dwarf elliptical experiencing a 5 mag- 
nitude brightening. Section 2 present the sample, observations 
and the methodology we followed. Kinematical and dynamical 
results are in section 3 and 4. We discuss our results in sec- 
tion 5 and a conclusion is given in section 6. In the following, 
we assumed a A-CDM cosmology with //q - 70, O,,, = 0.3, 
Ao = 0.7 and - -0.55. 

2. Sample, observations & methodology 

To select compact galaxies, we used Hubble Space Telescope 
(HST) images in the F8 14W filter (WFPC2, 0. 1 arcsec/pix; and 
ACS, 0.05 arcsec/pix) in the CFRS (3hr and 22hr fields) and the 
HDFS field. For one galaxy, however, we used ground-based 
images obtained by the CFRS team at the CFHT (0.207 arc- 
sec/pix, see Hammer et al. 2001). 

We selected 21 LCGs with 0.4 < z < 0.75, following the 
procedure detailed in Hammer et al. (2001), using the light con- 
centration parameter 5 as a compactness criterion which mea- 
sures the difference between the luminosities within the 5 and 
15 kpc radii (see Table [Q. The condition SI < 0.73 allows to 
select galaxies with rkaif < 5/igQ kpc in a homogeneous way. 
This value corresponds to ~ 4.34 kpc in a A-CDM cosmol- 

' They defined a compact galaxy as ri,aif < 0.5 arcsec and 56;8i4 < 
22.2 mag.arcsec"^. 



ogy at z ~ 0.55. To achieve a homogeneous selection, ACS 
images were degraded and re-sampled to the WFPC2 quality 
using a point spread function (PSF) generated with the Tiny 
Tim software. This size criterion ensures the selection of rela- 
tively small galaxies (see Ravindranath et al. 2004), although 
not necessarily as compact as those selected by some other au- 
thors (e.g. Guzman et al. 1997). Half light radii were then de- 
rived by interpolating luminosities enclosed within concentric 
ellipses using the IRAF polyphot task (see Hammer et al. 2001 
for a complete description of the procedure). Inclinations were 
estimated using Sextractor (Bertin et al. 1996) and the ellipse 
task of IRAF: we find a mean absolute difference of ~ 2 de- 
grees between these two methods. Independent measurements 
were also done by eye and gave similar results as Sextractor to 
within ~ 4 degrees. In the following, we will use the estimates 
from Sextractor and assume an error of +4 degrees. 

As part of the Guaranteed Time Observations (GTO) of 
the Paris Observatory, we obtained observations with the 
FLAMES/GIRAFFE instrument of the 2 1 compact galaxies us- 
ing the deployable integral field units (IFUs), each covering 
an area of 3 by 2 arcsec^, at 0.52 arcsec/pixel. The complete 
description of the GTO sample is given in Paper I (Flores et 
al. 2006). Briefly, we used LR04 and LR05 setups targeting 
the [Oil] doublet (R~ 10000), integration times ranging from 
8 to 13 hours and the seeing was typically ~ 0.6 arcsec dur- 
ing all the observations. Data cubes were reduced using the 
GIRBLDRS vl.l2 package (Blecha et al. 2000), including nar- 
row flat-fielding. Sky was carefully subtracted with our own 
IDL procedures. 

Among these 21 compact galaxies, we selected 17 galaxies 
for which at least 4 pixels had [Oil] doublet reaching a spectral 
signal to noise ratio (SNR) of 4 (see definition in Paper I). In 
the following we focus only on these 17 remaining LCGs (see 
Table^. At first sight, it might seem too challenging to derive 
velocity fields of compact galaxies (r/^;/ ~ 0.5 arcsec), using 
the GIRAFFE/IFU with 0.52 arcsec microlenses. To assess this 
we examine how many pixels that reach an integrated SNR of 
3. Among the sample of LCGs, we found a median value of 1 1 
pixels (~ 3 arcsec^) compared to 16 (~ 4.3 arcsec^) for a sam- 
ple of 8 spirals of Paper I (Flores et al. 2006). On average our 
LCGs are thus ~ 30% less extended than spirals. However, the 
average filling factor of the IFU (20 pixels) for the sample of 17 
LCGs is ~ 55 %, which is sufficient to explore the kinematics 
of these galaxies. 

We processed the spectra with a Savitzky-Golay filtering, 
which has the advantage over the widely used box smoothing 
that it conserves the first moment of spectral lines (Press et al. 
1989). We identified the [Oil] doublet by visual inspection and 
retained the spectra that reached a sufficient spectral SNR of 
3. We then fitted a double Gaussian with the following con- 
straints (where the subscripts denote the two components of 
the fit): /l2-/li=2.783A (in rest frame wavelength) and o-]-cr2. 
The line ratio was allowed to vary freely except when the fit 
failed: in these cases we forced the line ratio to a value of 1 .4 
which was the median value observed in our integrated spectra. 
This occurs for pixels with rather low SNR and affects ~ 10 % 
of the measured pixels and will thus not significantly affect our 
results. In all such cases, we checked by eye that the derived fit 
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was acceptable and took this into account during the classifica- 
tion (see below). The complete procedure is described in Paper 
I. 

Figure [T] shows velocity and velocity dispersion (cr) maps 
for the 17 LCGs. The cr -maps were coiTected for the instru- 
mental broadening using sky lines. To make the interpretation 
easier, velocity fields and cr-maps are presented after a simple 
5x5 linear interpolation. 

3. Kinematics of LCGs 

Following Paper I (Flores et al. 2006), we define 3 kinematical 
classes: rotating disks (velocity field showing rotation and cr- 
map showing a peak near the center), perturbed rotations (ve- 
locity field showing rotation but cr-map without peak or with 
a peak offset from the center) and complex kinematics, (see 
Figure[2and comments on individual objects). This classifica- 
tion relies on the fact that observations with a low spatial res- 
olution integral field spectrograph, such as GIRAFFE, should 
reveal line widths dominated by the integration of larger scale 
motions and not by intrinsic random motions: in the case of 
a rotating disk, cr should show a peak near the galaxy center, 
where the gradient of the rotation curve is the steepest (e.g. 
Van Zee & Bryant 1999). This classification has been checked 
through numerical simulations which are described in detail in 
Paper I. In these simulations, we assumed that all the observed 
galaxies are indeed rotating disks (taken as a standard model) 
and that all the observed large scale motions in the velocity 
fields correspond to rotations. In other words, we have tried 
to force each system to appear as rotational disks and then es- 
timate the discrepancy between the observed system and the 
adopted standard rotation (see Paper I). Below we comment on 
the velocity fields of individual objects. 

HDFS4170 

Its kinematics is classified as a rotating disk: a rotation is 
seen in the velocity field and the cr-map has a clear peak in 
the center. 

HDFS5190 

Its kinematics is classified as a rotating disk: a clear 
rotation is seen in the velocity field and the cr-map shows a 
peak near the center. 

CFRS03.0619 

Its kinematics is classified as a rotating disk: the velocity 
field shows a clear rotation and the cr-map has a peak 
slightly offset from the center of the galaxy. Due to the low 
spatial resolution of GIRAFFE, we choose to classify this 
galaxy as rotating disk rather than perturbed rotation. 

CFRS03.1032 

Its kinematics was first classified as complex because 
the velocity field looked perturbed and the cr-map has a 
peak near the edge of the galaxy. Although our simulation 
cannot reproduce the amplitude of the cr peak, this peak 
is however located in the same pixel as the one seen in 
the observed cr-map. HST/ACS imaging reveals a very 



compact structure, completely dominated by the center. 
This galaxy is the most compact in the sample. We 
retrieved a FORS slit spectrum from which we found a 
Log([OIII]/H^) ratio of 1.8, which characterize a Seyfert 
galaxy spectrum: the high cr value of this galaxy might 
reflects an AGN activity. We choose to re-classify its 
kinematics as a perturbed rotation, but we notice that it 
somewhat escapes our classification scheme. 

CFRS22.0619 

Its kinematics is classified as a perturbed rotation: a clear 
rotation is seen in the velocity field but the cr-map has a 
peak at the edge of the galaxy. Note that this galaxy is seen 
nearly edge-on. 

CFRS03.1349 

Its kinematics was first classified as a rotating disk because 
a clear rotation is seen in the velocity field and a peak is 
seen near the center of the cr-map. However, the simulation 
(see Paper I) cannot reproduce the location of the cr peak, 
although the secondary peak (at the bottom-left of the 
maximal peak) is reproduced. Zheng et al. (2004) classified 
this galaxy as Sab with a compact bulge, relatively blue 
compared to the disk. An interacting companion is 20 kpc 
away at the same redshift. Interestingly, this galaxy shows 
a companion galaxy at ~1 kpc, and the distorsion in the 
cr map is oriented towards this companion. We suspect 
that an interaction (gas falling) is responsible for both the 
star formation activity and distorsion of the kinematics 
(following the scenario proposed by Barton & Van Zee 
2001). We then choose to re-classify its kinematics as a 
perturbed rotation. 

CFRS22.1064 

Its kinematics is classified as perturbed: the velocity field 
shows rotation but the cr-map is very perturbed. 

HDFS5150 

Its kinematics is classified as a perturbed rotation: rotation 
is seen in the velocity field but the cr-map has a peak at the 
edge of the galaxy. The morphology looks quite irregular. 

CFRS03.0508 

The velocity field clearly shows an apparent rotation and 
the cr-map has a well-defined peak in the center. Note 
however that the dynamical axis seems almost orthogonal 
to the photometric axis of the brightest component: this 
could be a signature of outflows (Veilleux et al. 2005; 
Bosma, private communication). Hence its kinematics 
has been classified complex. From morphological studies, 
Zheng et al. (2005) classified this galaxy as the relics of an 
interaction or merger, with a relatively blue color over the 
whole galaxy. 

CFRS03.0645 

Also this galaxy shows rotation that is orthogonal to the 
photometric major axis, which may indicate an outflow 
and the cr-map has a peak at the edge of the galaxy. We 
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ID 


z 


1° 

AB 


K 


K 


SI 


half 




HDFS4170 


0.4602 


20.79 


-20.43 


-22.60 


0.56 


3.57 


51 


HDFS5190 


0.6952 


21.31 


-21.25 


-21.92 


0.60 


4.07 


59 


CFRS03.0619 


0.4854 


20.80 


-20.67 


-21.93 


0.63 


3.87 


27 


CFRS03.1032 


0.6180 


20.49 


-21.18 


-22.64 


0.51 


1.79 


37 


CFRS22.0619 


0.4676 


21.55 


-19.33 


-19.33 


0.70 


4.31 


68 


CFRS03.1349 


0.6155 


20.87 


-21.19 


-22.92 


0.56 


3.84 


48 


CFRS22.1064 


0.5383 


22.08 


-19.87 


-21.64 


0.24 


2.36 


49 


HDFS5150 


0.6956 


22.36 


-20.20 


-21.03 


0.72 


3.39 


42 


CFRS03.0508 


0.4642 


21.92 


-19.61 


-20.35 


0.46 


3.32 


38 


CFRS03.0645 


0.5275 


21.36 


-20.30 


-21.34 


0.71 


4.57 


45 


CFRS22.0919 


0.4738 


21.77 


-19.99 


-19.53 


0.41 


2.52 


61 


CFRS22.0975 


0.4211 


20.21 


-20.40 


-22.53 


0.59 


3.82 


50 


CFRS03.0523 


0.6508 


21.31 


-20.67 


-21.55 


0.48 


3.57 


41 


HDFS4130 


0.4054 


20.09 


-20.90 


-22.13 


0.62 


4.04 


36 


HDFS4090 


0.5162 


22.15 


-19.71 


-19.83 


0.35 


1.54 


45 


HDFS5140 


0.5649 


22.38 


-19.76 


-20.46 


0.36 


2.57 


50 


HDFS5030 


0.5821 


20.40 


-21.74 


-22.68 


0.66 


4.19 


25 



" isophotal magnitudes, 
from Hammer et al. 2005. 
in kpc. 

inclination, in deg. 

Table 1. Main properties of the sample of LCGs: galaxies names, redshifts, isophotal I magnitudes, absolute B and K magnitudes, 
light concentration, half light radii, inclination. 



classified the kinematics as complex. Zheng et al. (2005) 
found a relatively blue color all over the galaxy and 
classified it as a probable merger. 

CFRS22.0919 

Its kinematics is classified as complex: the velocity field is 
perturbed and the cr-map does not show any peak. Note the 
tails, characteristic of interacting systems, seen in the HST 
image. 

CFRS22.0975 

Its kinematics is classified as complex: the velocity field 
is perturbed and the cr-map shows a peak but not at the 
galaxy center HST imaging reveals 3 distinct components. 
This system is probably just preceding a merger and the 
high velocity gradient in Figure [2 can probably not be 
attributed to rotation. 

CFRS03.0523 

Its kinematics is classified as complex: the velocity field 
and the cr-map do not show any kind of structure expected 
from rotation. HST imaging shows a tidal tail, probably 
indicating ongoing interactions or gas accretion. This 
galaxy has a central region bluer than the outer region 
(Zheng et al. 2004). The cr maxima on the outer edges of 
the galaxy correspond to relatively blue regions (see Zheng 
et al. 2004). 

HDFS4130 

Its kinematics is classified as complex: the velocity field 
shows rotation but the cr-map does not have any peak. This 



galaxy is asymmetric with (spiral ?) arms visible on one 
side only (maybe a distant version of a "tadpole" galaxy). 

HDFS4090 

Its kinematics is classified as complex: the velocity field 
is perturbed and the cr-map does not have any peak. The 
morphology looks like a "peanut". 

HDFS5140 

Its kinematics is classified as complex: the velocity field 
is very perturbed, although the cr-map has a peak near the 
center 

HDFS5030 

Its kinematics is classified as complex: both the velocity 
field and the cr-map appear perturbed. 

In the sample, 53% of the LCGs (9 galaxies) present com- 
plex kinematics, which are very different from classical nearby 
spiral galaxies (e.g. Garrido et al. 2002; 2003 & 2004) or early- 
type galaxies (Emsellem 2004), indicating that these LCGs are 
unrelaxed systems. As Figure [2 shows, there is no correlation 
between these dynamical classes and their SNR. Even discard- 
ing the lowest SNR pixels would not change the nature of the 
kinematics and hence the kinematical complexity of a galaxy 
cannot be attributed to a lower SNR. 

Five galaxies, or 29%, have perturbed kinematics and 
the remaining three galaxies (18%) look like normal spiral 
galaxies. Interestingly, among these last eight galaxies (per- 
turbed rotation and rotating disks), all but three (CFRS03.0619, 
HDFS5150 and HDFS5190) have possible nearby companions 
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which might indicate interactions responsible for their com- 
pactness (see Barton & van Zee 2001). 

Before concluding on the dynamical nature of LCGs, we 
have to investigate whether such perturbed/complex velocity 
fields could be artificial features caused by the low spatial sam- 
pling of GIRAFFE. In Flores et al. (2004 & 2006), we illus- 
trated the ability of GIRAFFE to recover regular velocity fields 
of distant rotating disks. The question is then to demonstrate 
the ability of GIRAFFE to recover also more complex velocity 
fields in distant galaxies. To tackle this issue, we used Perot- 
Fabry observations by Ostlin et al. (1999) of a local LCG. We 
redshifted the ESO400-G43 data cube (see Ostlin et al. 1999 
& 2001) to z=0.6 and simulated the effects of a 0.52 arcsec 
sampling and a 0.6 arcsec seeing. The result is illustrated in 
Figure 12 GIRAFFE smoothes the velocity fields (and thereby 
underestimates the velocity gradient) but no artificial features 
are created. Moreover, the simulated velocity field shows some 
similarities with those of the LCGs in Fig. 1. 




Fig. 2. Simulation of the velocity field of a LCG at z=0.6. 
Left: Redshifted velocity field of ESO400-G43 (see Ostlin et 
al. 1999) with a ~ 0.1 arcsec pixel (1/5 of a GIRAFFE pixel). 
Right: velocity field of the same galaxy as could have been seen 
by GIRAFFE (0.6 seeing and 0.52 arcsec/pixel) with a simple 
5x5 interpolation and assuming the same spatial extension. The 
resolution is dramatically reduced as well as the velocity gra- 
dient (by a factor 4.7), but no artificial features are created. 

We can then conclude that LCGs are mainly (53 to 82% 
as judged from our rather small number statistics) unrelaxed 
systems, at least when ionized gas motions are considered. 
Such kind of complex kinematics have also been observed by 
Swinbank et al. (2005) who obtained the velocity field of a 
merger remnant at z ~ 0.1, similar to the complex velocity 
fields of our LCGs. However, this conclusion relies on the kine- 
matics of the gas only: Ostlin et al. (2004) showed that in one 
local LCG, gas and stars seem dynamically decoupled, and that 
stars seem more relaxed than the gas. In the following section, 
we attempt to investigate the nature of the LCG kinematics, 
namely if they are mostly dominated by rotation or by disper- 
sion. The fact that LCGs are not dynamically relaxed systems 
is a severe limitation for this exercise. 

4. Test of the dynamics of LCGs 

To investigate energy balances of LCGs, we assume in this 
section that LCGs are systems almost at equilibrium. In this 
case, one could imagine two possibilities: (I) LCGs are mainly 
supported by velocity dispersion, or (2) LCGs are mainly sup- 
ported by rotation. Both options require to assume that LCGs 



are dynamically not too far from equilibrium and we would 
like to point out again that this assumption represent a severe 
limit in setting up energy balances of LCGs. Moreover, the last 
possibility (that LCGs are dominated by rotation) requires in 
addition to identify large scale motions to rotation, which is 
even more speculative and wiU thus be detailed in appendix. 

4.1. Can LCGs be supported by velocity dispersion? 

If LCGs were systems at equilibrium and supported by velocity 
dispersion, then the intensity-weighted mean cr and the one de- 
rived from their integrated spectra cr,„,y should be in agreement 
(e.g. Bershady et al. 2004; Ostlin et al. 2001). 

We constructed integrated spectra by the direct summation 
of all the spatial elements in the whole IFU. Integrated [Oil] 
lines were fitted following the same way as for IFUs pixels, but 
due to the influence of larger scale motions (rotation for spirals) 
that widen integrated lines, about 50% of the integrated spectra 
were impossible to fit coiTectly by a double Gaussian. We then 
summed independently the fits of each line of the [Oil] doublet 
over the whole IFU, fitted both by single Gaussian (corrected 
from instrumental dispersion) and estimated the integrated ve- 
locity dispersion taking the mean of the two velocity disper- 
sions derived independently. We checked that in most cases, 
these two velocity dispersions were very similar and that when 
integrated line fit was possible, both methods gave similar re- 
sults. 

Both measurements are presented in Table |3 and Figure 
13 We estimated the error on the sigma measurement to 
be 10% (median, see Paper I). We thus adopt a 30% 
relative treshold (3-sigma treshold) between cr and o",-,,,;, 
to identify the galaxies which could potentialy be sup- 
ported by velocity dispersion. Nine galaxies (roughly 50%) 
have such a relative difference between cr and cr,„,^: 
CFRS03.0523, CFRS03.1032, CFRS22.0919, HDFS4090, 
HDFS5030, HDFS5150, CFRS03.0645 CFRS22.1064 and 
HDFS4130 (see Tableland Figure|3. Note that among these 
galaxies, CFRS03.1032 is a very peculiar case because of both 
its very high central cr and its AGN activity (see individual 
comments). 

4.2. Dynamical support of LCGs 

The other possibily for LCGs to be systems at equilibrium is to 
be supported by rotation. This would imply that the large scale 
motions seen in the velocity fields of LCGs are due to rotation, 
althgough these motions are not completely relaxed. However, 
we already know that this is not true for at least the objects 
whose dynamical axis are not aligned with the optical one, ie 
for objects suspected of outflows (see individual comments). 
In the appendix, we nevertheless naively assumed that the large 
scale motions are precisely due to rotation. As we assumed that 
LCGs are at equilibirum, we can then set up energy balance us- 
ing classical relations linking their kinematics (velocity or dis- 
persion) to their mass, to see if these balances are effectively 
dominated by the large scale motions interprated as rotation. 
Under such an assumption, one finds that 70% of LCGs have an 
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ID 


class 


V,-o, 


O" 


C intg 


M, 


HDFS4170 


RD 


173 


75 


111 


10.83 


HDFS5190 


RD 


168 


64 


126 


10.51 


CFRS03.0619 


RD 


155 


57 


88 


10.50 


CFRS03.1032 


PR* 


139 


162 


168 


10.87 


CFRS22.0619 


PR 


67 


36 


64 


9.50 


CFRS03.1349 


PR* 


235 


70 


183 


10.91 


CFRS22.1064 


PR 


124 


107 


124 


10.39 


HDFS5150 


PR 


77 


69 


78 


10.14 


CFRS03.0508 


CK 


78 


37 


54 


9.87 


CFRS03.0645 


CK 


125 


58 


67 


10.27 


CFRS22.0919 


CK 


45 


46 


51 


9.54 


CFRS22.0975 


CK 


409 


82 


182 


10.82 


CFRS03.0523 


CK 


97 


156 


148 


10.35 


HDFS4130 


CK 


153 


76 


97 


10.63 


HDFS4090 


CK 


20 


51 


49 


9.67 


HDFS5140 


CK 


224 


69 


133 


9.91 


HDFS5030 


CK 


74 


59 


65 


10.81 



Table 2. Dynamical properties of the sample of LCGs. The column entries are (from left to right): id, dynamical class 
(RD=Rotating Disk, PR=Perturbed Rotation, CK=Complex Kinematics), maximal rotational velocity (corrected from inclina- 
tion but not from GIRAFFE spatial sampling; The objects for which the kinematical classification is uncertain are indicated by a 
star, see individual comments), intensity-weighted velocity dispersion calculated from cr-map, velocity dispersion derived from 
integrated spectra and stellar masses derived from photometry (in solar mass). 
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Fig. 3. Velocity dispersion o",-,,,;, derived from integrated spec- 
tra vs intensity weighted mean velocity dispersion cr derived 
from cr-maps. Encircled blue full dots are LCGs classified as 
rotating disks, encircled green squares are LCGs classified as 
perturbed and encircled red triangles are complex LCGs. Blue 
full dots are rotating disks from Paper I (Flores et al. 2006), 
added for comparison. The dash line represents the limit where 
the relative difference between cr and cr is 30%. 



energy balance effectively consistent with rotation. However, 



this question can only be properly addressed with the knowl- 
edge of the kinematics of the stars (see e.g. Ostlin et al. 2004). 

To summarize, we find that about 50% of LCGs could be 
supported by velocity dispersion, assuming they are not too far 
from equilibrium. Given that 18% of LCGs are classified as 
RD, it remains one third of LCGs for which we cannot exclude 
that a rotational support could play a role in their dynamical 
state. Both large scale and random motions seem to play an 
important role in about 40% of LCGs. Such a mix is perfectly 
compatible with mergers. 

5. Stellar vs. dynamical masses 

We derive stellar masses M» following Hammer et al. (2005) 
(see Table 2). Figure |4] compares stellar masses with total 
dynamical pseudo-equivalent masses derived from integrated 
spectra, which is the most reliable estimate of the total dynam- 
ical mass we can use (see Appendix). Following Conselice et 
al. (2005) who carried an extensive study of uncertainties and 
systematics on stellar masses, we estimate our error on stellar 
masses to be 0.5 dex. Among the sample of spirals of Paper I 
(Flores et al. 2006), a mean ratio of dynamical to stellar masses 
of ~ 6 (median ~ 7) was found, which is roughly in agreement 
with the value of Conselice et al. (2005) who found a mean ra- 
tio of ~ 5 for a sample of z ~0.5 disk galaxies. The difference 
between these two samples of spirals is most probably related 
to the difference of observational strategy (slit spectroscopy vs 
Integral Field Spectroscopy), because the maximal rotational 
velocity can be largely underestimated when measured by slit 
spectroscopy, since only a part of the whole kinematics is sam- 
pled (see Discussion). On the other hand, we find for our sam- 
ple of LCGs a mean ratio between dynamical and stellar 
masses of ~ 2 (median ~ 1.6). If the majority LCGs are merg- 
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ing systems, as suggested by their kinematics, their spectra are 
then Hkely dominated by unrelaxed motions, which could lead 
to underestimate their total dynamical Ai^. . 




log ME /Mg 
iiiig 

Fig. 4. Stellar masses vs dynamical Ai^ . Encircled blue full 
dots are LCGs classified as rotating disks, encircled green 
squares are LCGs classified as perturbed and encircled red tri- 
angles are complex LCGs. Blue full dots are rotating disks from 
Paper I (Flores et al. 2006), added for comparison. For conve- 
nience, only the median error bar has been plotted. 



6. Discussion & conclusion 

We have presented GIRAFFE-IFU observations and derived 
velocity fields and cr maps for 17 LCGs. Our main result is 
that 18% have velocity fields characteristic of rotating disks 
kinematics, the rest having perturbed or complex kinematics. 
This result is unlikely to be affected if we account for the four 
objects that were discarded because of a too low SNR. Among 
them, two show (in HST images) tails characteristics of ongo- 
ing interactions. Assuming that they are not too far from equi- 
librium, we find that about half of LCGs could be supported 
by velocity dispersion. The remaining 32% of LCGs seems to 
be dominated by large scale motions. An important fraction of 
LCGs, between 32 and 82%, are thus probable mergers. 

In this sample of 17 LCGs, four have their velocity gradient 
nearly perpendicular to their main optical axis (CFRS03.0508, 
CFRS03.0645, CFRS22.0975 and HDFS5140). It suggests that 
in these galaxies, gas is tracing feedback processes such as 
outflows, rather than gravitational dynamics: in this case, the 
dynamical axis should appear off axis relative to the optical 
main axis (see Veilleux et al. 2005 and references therein). All 
of them except CFRS22.0975 have stellar masses lower than 
3.1O'"M0 which is the upper threshold for supernova feed- 



back to drive efficient outflows (Dekel & Birnboim 2005). In 
CFRS22.0975, the gas velocity field probably traces the rela- 
tive motion of the merger progenitors (see its morphology in 
Figure 1). For CFRS03.0645, we retrieved FORS slit spectra 
(Hammer et al. 2001) to compare systemic velocities of emis- 
sion (H/3, Hy, [Oil] and [OIII]) and absorption Unes (Call, H 
and K) and found no significant shifts between them, which 
makes the outflow hypothesis uncertain for this galaxy. Note 
however that the slit has been aligned with the apparent (op- 
tical) main axis which could explain why we did not see any 
shift between emission and absorption lines. The electron den- 
sity map of CFRS03.0508 is presented in Paper III (Puech et 
al. 2006b) and support the outflow hypothesis. 

Bershady et al. (2004) used the STIS long slit spectro- 
graph onboard the HST to study the kinematics of 6 LCGs with 

Mb 21, rhaif ~ 2 /i^g kpc and cr ~65 km/s, and found that 

LCGs are supported by velocity dispersion. Unfortunately, only 
one galaxy of the present sample (CFRS220919) is enough 
compact to fulfill these selection criteria. Interestingly, it has 
being pointed out by Hammer et al. (2001) that this galaxy 
could be a possible progenitor of a dE. Here we find that this 
galaxy is a possible candidate for being supported by disper- 
sion. However, given the complexity of most LCGs kinemat- 
ics presented here, slit spectroscopy should be used with care 
for this kind of objects. A careful inspection of Figure ^ re- 
veals that with long slit spectroscopy, the real nature of many 
LCGs presented in this paper would have been misinterpreted. 
In the case of CFRS03.0645, a slit positioned along the main 
optical axis would have revealed a rather flat velocity gradi- 
ent and would have completely missed the rotation. The case 
of CFRS03.0508 is even more instructive: the same exercise 
would have revealed a flat velocity gradient and a clear cr peak, 
and this galaxy would then have been classified as supported 
by dispersion. 

Finally, how can we interpret the compactness of LCGs? 
Most of perturbed kinematics and rotating LCGs show possible 
companions with which they could be in interaction. This could 
explain their compactness as due to interactions and/or minor 
mergers, following Barton & van Zee (2001, see introduction). 
Another conjecture is the one of Hammer et al. (2005). They 
proposed a scenario where local massive spiral could form af- 
ter major mergers in three main phases. The sequence would 
start by a pre-merging phase during which the system would 
form a huge amount of stars and appear as a LIRG. The sec- 
ond phase would be the LCG phase, were all material falls 
onto the mass barycenter of the merging system which could 
enhanced the star formation activity in the center of these sys- 
tems, making them looking compact. This is consistent with 
Bergvall & Ostlin (2002) who found central intense starbursts 
superimposed on low surface brightness components in four 
local BCGs. Our results are consistent with this picture as we 
found that most LCGs are objects with complex kinematics as 
expected from major mergers. During the third phase of the 
scenario, a disk would grow thanks to material accreted from 
the IGM. Inflows/outflows are also predicted by this scenario, 
arising from feedback and gas falling back to form a new disk 
(see Robertson et al. 2005). 



8 



M. Puech, F. Hammer, H. Floras et al.: Are LCGs merger remnants ? 



How could we distinguish between minor and major merg- 
ers? In a minor merger, the disk cannot be destroyed and the 
kinematics of the remnant cannot appear too complex. We 
should then observe a galaxy still rotating along its main optical 
axis but a dispersion map peaking outside the center, where the 
smaller progenitor falls. This could correspond to the LCGs we 
classified perturbed. On the other hand, during a major merger 
the disk is completely destroyed or at least strongly perturbed. 
In such a case, we should then observe either a rotation signif- 
icantly misaligned with the optical axis combined with a non- 
centered dispersion map, or a complex kinematics without any 
obvious structure. We emphasize that given our spectral resolu- 
tion, reaching ~ 10000, the fact that the [Oil] doublet is not al- 
ways resolved reveals by itself the complexity of some of these 
galaxies. 

LCGs dominate the evolution of the UV luminosity from 
z =1 to 0. Their role can thus not be negligible in the process of 
formation and evolution of galaxies during the last 8 Gyr Our 
results highly support a hierarchical type picture where galax- 
ies form from smaller units. In this picture, LCGs seem to be 
a major event as proposed by Hammer et al. (2005). Although 
based on small numbers, our sample is nevertheless represen- 
tative of the galaxy population at 0.4 < z < 0.75 (see Paper 
I). A larger sample is under construction as part of the ESO 
Large Program IMAGES (P.I.: F. Hammer) and will be a deci- 
sive step towards the confirmation of the spiral rebuilding sce- 
nario. Recently, it has been shown that major merger remnants 
may not necessarily be ellipticals but also spirals, depending 
on the gas abundance (Springel and Hernquist 2005) and/or 
the Black Hole feedback (Robertson et al. 2005). Comparisons 
with theoretical simulations will bring a crucial test of the spi- 
ral rebuilding scenario and on the nature of LCGs. 
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Appendix A: Energy balance of LCGs 

This appendix is devoted to the energy balance of LCGs. It is 
explicitly assumed that LCGs are systems at equilibrium and 
supported by rotation. We will thus assume in this appendix 
that the large scale motions seen in the velocity fields are asso- 
ciated with rotation, even if the true origin of these large scale 
motions is largely uncertain, except for some suspected out- 
flows and those clearly identified as RD from their kinematics 
(see individual comments and section 4.2). We will neverthe- 
less naively assume that their origin is rotation and see if any 
contradiction arises. To test this hypothesis, we set up an energy 
balance in the sample of LCGs. In this balance, we will take 
into account the contributions from large scale ordered motions 
(interpreted as rotation) and from random motions. Energies 
are estimated in what we call "pseudo equivalent masses" ( M^, 
i.e. in mass units). It is important to emphasize that these Ai^ 
can be interpreted as real masses only for galaxies in equilib- 
rium. 
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A.1. Contribution from rotation 



To estimate the mass supported by rotation, we assumed that 
the maximal rotational velocity is equal to half the maximal 
gradient of the velocity field, corrected for inclination (see 
Table 1). Due to the low spatial sampling and large distances to 
our targets, we have already pointed out that the GIRAFFE/IFU 
observations will underestimate the maximal rotational veloc- 
ity (see Figure Using hydrodynamical simulations of an 
Sbc (Milky-Way like) galaxy by Cox et al. (2004), we simu- 
lated our GIRAFFE observations assuming median atmospher- 
ical conditions at ESO VLT (0.81 arcsec seeing at 500 nm and 
an outer scale for the turbulence of 24 m, see Tokovinin 2002). 
We scaled this template galaxy (Vroi - 160 km/s and inc = 53 
degrees) to fill boxes of length ranging from 0.75 to 6 arcsec to 
mimic distant galaxies, and then compared the kinematics seen 
by GIRAFFE with the original simulation (see Figure lATl . We 
found that for spiral galaxies with sizes between 2 and 3 arcsec, 
GIRAFFE is able to correctly recover the maximal rotational 
velocity, although it is underestimated by ~ 20 %. For more 
complex kinematics, the correction factor should be larger (see 
also Figure|2}: between 2 and 6. However, as we explicitly as- 
sume in this section that LCGs are rotating disks, we applied 
a constant factor of 20% whatever the dynamical class (spiral 
or perturbed/complex) of LCGs in order to compute homoge- 
neous estimates. 




12 3 4 



Galaxy size (arcsec) 

Fig. A.l. Correcting factor for the maximal rotational velocity 
vs galaxy size. Due to its low spatial sampling, the GIRAFFE 
IFU leads to underestimate Vroi by a factor 1 .2 + 0.04 for galax- 
ies sizes between 2 and 3 arcsec. Full line: tilted view (~ 53 
degrees). Dash hne: top view (~ 25 degrees) 



Rotational are then estimated from (Lequeux 1983): 



with f=0.6 for a disk with a flat rotation curve and R25 - 
'^f'haif, following Phillips et al. (1997). For a pure exponential 
disk with scale length r,/, r/,,,// = 1.67 rj. Taking R25 = 2rhaif 
in fact assumes that rhaif/rd = 1.6 since R25 = 3.2rd (Persic & 
Salucci 1988). With in Mq, Vrot in km/s and r/,a// in kpc, 
this becomes: 

Mf = 0.279 • 10V,,„,/y2„, 

Corrected rotational velocities and Ai^ are given in Table 
2. To estimate our error bars we took into account the uncer- 
tainties on rfiaif (~l/2 HST pixel which represents 0.34 kpc at 
z=0.6), Vrot (± 5 km/s, determined by repeating several times 
the fitting procedure), inclination (+4 degrees, see section 2) 
and the correction factor on the velocity (+0.04, see Figure 
IA.1> . Note however, that the real uncertainty on inclination is 
probably higher for unrelaxed systems, as its derivation usually 
relies on the assumption of a thin disk seen in projection (see 
section 2). From this we estimated an median error on Alf of 
~ 0. 16 dex which is mainly dominated by the uncertainty in the 
inclination. 

A.2. Contribution from velocity dispersion 

Following Ostlin et al. (2001), we estimated the supported 
by velocity dispersion using (M^ in Mq, cr in km/s and r/,;,// 
in kpc): 

M^ = l.l-10V,w/cr2 

where cr is the intensity weighted mean of the cr -map. Due to 
the low GIRAFFE spatial sampling, cr will tend to be overes- 
timated, since the coarse pixel size integrates large scale mo- 
tions (rotation for a spiral galaxy). We used the same hydrody- 
namical simulation to estimate the increase in cr due to large 
scale motions and found ~ 50 km/s {Vrot - 160 km/s and 
inc - 53 degrees). We used another simulation to check that 
the effect approximatively scales with 1 / sin( /nc) (within an er- 
ror on cr of 5 km/s) and assumed the same scaling with Vrot- 
We then corrected the GIRAFFE intensity weighted mean cr 
using this recipe and from this estimated the mass supported 
by velocity dispersion (see Table lA.ll and Figure IA.2> . For 
pure rotating disks, a correction roughly equal to the measure 
is expected, which means that intrinsic velocity dispersion is 
negligible. Galaxies classified as rotating disk and which fall 
far from the cr - crdiskmodei region could be galaxies with 
a significant bulge (e.g. CFRS03.9003, see also Puech et al. 
2006b) or show perturbation in their cr-map (CFRS03.0619, 
CFRS22.0504 and HDFS4020). Interestingly, these 4 galaxies 
(plus CDFS03.1032, but see individual comments) are among 
those which are the nearest of the cr;„,^ - cr line on Figure 
|3] which supports the idea that their dynamical support could 
include a substantial contribution coming from dispersion. All 
points with a correction larger than the measured mean cr are 
spiral galaxies except one (CFRS22.0975) whose velocity gra- 
dient is likely due to relative motions between merging compo- 
nents (see its morphology in Figure rather than rotation. In 
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ID 


class 




Mf 


Ml 






HDFS4170 


RD 


56 


10.63 


10.10 


10.74 


10.69 


HDFS5190 


RD 


27 


10.66 


9.50 


10.69 


10.85 


CFRS03.0619 


RD 


55 


10.57 


10.10 


10.70 


10.52 


CFRS03.1032 


PR* 


161 


10.14 


10.71 


10.81 


10.75 


CFRS22.0619 


PR 


24 


9.89 


9.43 


10.02 


10.29 


CFRS03.1349 


PR* 


34 


10.93 


9.70 


10.94 


11.15 


CFRS22.1064 


PR 


101 


10.17 


10.43 


10.62 


10.60 


HDFS5150 


PR 


67 


9.90 


10.23 


10.40 


10.35 


CFRS03.0508 


CK 


34 


9.91 


9.63 


10.09 


10.03 


CFRS03.0645 


CK 


49 


10.46 


10.09 


10.60 


10.35 


CFRS22.0919 


CK 


43 


9.30 


9.70 


9.85 


9.85 


CFRS22.0975 


CK 


10 


11.41 


8.62 


11.41 


11.15 


CFRS03.0523 


CK 


155 


10.13 


10.97 


11.03 


10.93 


HDFS4130 


CK 


72 


10.58 


10.36 


10.79 


10.62 


HDFS4090 


CK 


51 


8.38 


9.64 


9.66 


9.61 


HDFS5140 


CK 


32 


10.71 


9.45 


10.74 


10.70 


HDFS5030 


CK 


59 


9.96 


10.21 


10.40 


10.29 



Table A.l. Dynamical properties of the sample of LCGs. The column entries are (from left to right): id, dynamical class 
(RD=Rotating Disk, PR=Perturbed Rotation, CK=Complex Kinematics), intensity-weighted velocity dispersion corrected from 
GIRAFFE sampling effect (see text), supported by rotation (corrected from GIRAFFE sampling effect), supported by 
random motions, total dynamical , and dynamical estimated by integrated velocity dispersion 



these few cases, we fixed the corrected value to 10 km/s which 
corresponds to the minimal expected line width due to intrinsic 
turbulent motions in spiral galaxies (Rozas et al. 1998; Van Zee 
& Bryant 1999). We estimate a median error on of ~ 0.08 
dex. 
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Fig. A.l. Correction term (Tdiskmodei vs intensity weighted mean 
cr. Encircled blue full dots are LCGs classified as rotating disks, 
encircled green squares are LCGs classified as perturbed and 
encircled red triangles are complex LCGs. Blue full dots are 
rotating disks from Paper I (Flores et al. 2006), added for com- 
paiison. 



A.3. Total dynamical 

To estimate total dynamical pseudo-equivalent masses Al^^,^, 
we simply added Aff and M^. Most authors derive cr on inte- 
grated spectra (global velocity dispersion) and use the relation 
of the previous section to estimate the whole dynamical mass 
Ml^^^^ (Guzman et al. 1996, Phillips et al. 1997). Under as- 
sumptions concerning the anisotropy of the kinetic energy ten- 
sor and the geometry of the system, this relation can be used 
to estimate total dynamical masses of rotating and flattened 
spheroids (Bender et al. 1992). This approach is often used in 
studies using slit spectroscopy (Guzman et al. 1997; Guzman 
et al. 2001; Hammer et al. 2001) but the validity of this relation 
for systems dominated by rotation is uncertain. In the follow- 
ing, we also estimated A1^,_^^ which we compared with M^^^. 

Figure IA.3I shows the comparison between M^^,^ and 
. We find a correlation between the two estimates, which 
seem to validate our estimates and is consistent with the fact 
that LCGs could be systems not too far from equilibrium. Spiral 
galaxies are almost equally distributed on both sides of the line 
where - M^, which is likely due to the fact that 

does not correct explicitly for inclination effects. 

In Figure IA.4I we compare My with M^^,^^ to test if the 
LCGs are dominated by rotation. For comparison, we also 
plot a sample of 8 rotating disks taken from Paper I (Flores 
et al. 2006). Clearly, the whole spiral sample falls in the ro- 
tation dominated area. At least 70% of the LCGs (12 galax- 
ies) seems to be dominated by rotation whereas the remaining 
30% (5 galaxies) seems to be dominated by velocity dispersion. 
Note that these five galaxies were already identified as poten- 
tially supported by dispersion (see section 4. 1): CFRS22.0919, 
CFRS03.0523, CFRS03.1032 HDF4090 and CFRS22.1064. 
Interestingly, one of these galaxies (CFRS22.0919) was iden- 
tified by Hammer et al. (2001) as a potential dwarf progenitor 
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Fig.A.3. Dynamical derived from integrated spectra vs 
dynamical Ai^ derived from 3D spectroscopy. Encircled blue 
full dots are LCGs classified as rotating disks, encircled green 
squares are LCGs classified as perturbed and encircled red tri- 
angles are complex LCGs. Blue full dots are rotating disks from 
Paper I (Flores et al. 2006), added for comparison. For conve- 
nience, only the median error bar has been plotted. 

as discussed by Guzman et al. (1997). Because we are likely 
underestimating their rotational velocities, we cannot exclude 
that all LCGs have energy balance dominated by rotation. 

However, we recall that this conclusion relies on the as- 
sumptions that large scale motions in the velocity fields are due 
to rotation, and that we know that for some systems (e.g. the 
suspected outflows) this is far from being true. We nevertheless 
cannot exclude that a possible relation could link most of the 
LCGs with rotating disks. 
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Fig.A.4. Rotational Ai^ vs total dynamical M^. Encircled 
symbols are the sample of LCGs. Encircled blue full dots are 
LCGs classified as rotating disks, encircled green squares are 
LCGs classified as perturbed and encircled red triangles are 
complex LCGs. Blue full dots are rotating disks from Paper 
1 (Flores et al. 2006), added for comparison. The dash line rep- 
resent the limit where M^-Ai^^.J2. For convenience, only the 
median error bar has been plotted. Notice that the high Alf 
of CFRS22.0975 (the point at the top-right) could not be at- 
tributed to a normal rotation since it is an obvious merger. 
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Fig.l. SEE ATTACHED JPEG FIGURE Kinematics of 17 LCGs. From left to right: I band HST imaging (0.1 arcsec/pix, 
FoV=3x2 arcsec^) with the IFU bundle superimposed; velocity field (5x5 interpolation, ~ 0.1 arcsec/pix), cr-map (5x5 interpo- 
lation, ~ 0.1 arcsec/pix) and SNR map (mean spectral SNR, see Paper I). From up to down: LCGs classified as rotating disks 
(3 first), perturbed (5 next) and complex (9 last). The 3 first galaxies classified as complex (CFRS03.0508, CFRS03.0645 and 
HDFS5 140) are those suspected to be dominated by outflows. 



Fig. 1. SEE ATTACHED JPEG FIGURE continued. 



This figure "lcg_final_partl.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/ astro-ph/0603 5 62 v 1 



This figure "lcg_final_part2.jpg" is available in "jpg" format from: 



http ://arXiv . org/ps/ astro-ph/0603 5 62 v 1 



